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1 Introduction 



The Bethe Ansatz solution of the open spin-| XXZ quantum spin chain with diagonal bound- 
ary terms has long been known P1I2]- However, the case of nondiagonal boundary terms [3] 
has resisted solution for many years (see, e.g., jl]). A Bethe Ansatz solution for the latter 
case has recently been proposed in [HI EJ (see also [Zj). In terms of the parameters introduced 
there, the Hamiltonian is given by 

1 JV_1 

H = 2 { S («+! + ^n+l + COSh V 



n=l 



+ sinh rj |^ coth a_ tanh f3-o{ + csch a_ sech /3_ ( cosh 9-<j\ + i sinh 9-<yf) 

— coth a + tanh (3 + a z N + csch a + sech f3 + ( cosh 9 + cr% + i sinh O+a^) }; (1-1) 

where a x , a y , a 2 are the usual Pauli matrices, r\ is the bulk anisotropy parameter, a±, /3± , 9± 
are boundary parameters, and N is the number of spins. An unusual feature of this Bethe 
Ansatz solution is that the boundary parameters must satisfy the linear relation 

a_ + /?_ + a+ + /3+ = ±(0_ - 9+) + rjk , (1.2) 

where k is an even integer if N is odd, and is an odd integer if N is even. 1 The energy 
eigenvalues are given by 

M x x 

E = sinh 2 77 > — -H — sinh rj (coth ct_ + tanh (3- + coth a + + tanh (3+) 

sinh Uj sinh(«j + 77) 2 



+ ~(JV- 1) cosh 77, (1.3) 



M, (1.4) 



where the Bethe roots {uj} satisfy the Bethe Ansatz equations 

h(-Uj - 7]) Q(Uj - 7]) ' 

where h(u) is given by 2 

7 / x • ,2N ( . xSinh(2M + 2r7) 

= -smh (w + r/) , 

smn(z« + 77) 

x 4 sinh(u + a J) cosh(w + (3 J) sinh(u + a + ) cosh(w + f3 + ) , (1.5) 



1 An alternative solution was proposed in [H] which does not require any constraint among the boundary 
parameters. However, that solution holds only for 77 values corresponding to roots of unity, and the Bethe 
Ansatz equations are not of the conventional form. 

2 This expression for h(u) differs from (3.26) in j^j by the factors as the result of working here with 

a rescaled transfer matrix, as discussed further in Section l2~Tl 
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and Q(u) is given by 

M 

Q(u) = Y\ sinh(w — Uj) sinh(u + Uj + r/) , (1.6) 

i=i 

which satisfies Q(u) = Q(—u — if). The Bethe Ansatz equations assume a more symmetric 
form if expressed in terms of the shifted Bethe roots 

% = ^ + |- (1-7) 

Another unusual feature of this solution is that the number M of Bethe roots is fixed for 
given values of N and k (similar to the case of the XYZ chain), and is given by 

M=i(JV-l + fc), (1.8) 

where k is the integer appearing in (jl.2jl . 

Several important issues were left unresolved in t>j. In particular, the Bethe Ansatz 
solution was obtained for the specific values of the anisotropy parameter rj = m /2 , iir/A , . . . 
(for which q = e v equals certain roots of unity). The solution was conjectured to hold for 
generic values of rj, but little direct evidence was given. Also, because the value of the integer 
k in the constraint (jl.2j) is tied to the number of Bethe roots through (J1.8)) . it is evident 
that the requirement of completeness should restrict the value of k. However, the problem 
of determining those restrictions remained unsolved. 

The purpose of this article is to address these and related questions. From numerical 
studies of chains with sizes up to N = 7, we find significant evidence that the Bethe Ansatz 
solution indeed holds for generic values of rj. In particular, for generic values of both bulk and 
boundary parameters, the Bethe Ansatz solution yields the complete set of 2 N eigenvalues 
when k = N + 1 (i.e., M = N). While it is gratifying to obtain all the eigenvalues, this result 
is also disappointing, as it is impractical to satisfy the constraint (jl.2|) with k = N + 1 in 
the thermodynamic (N — > oo) limit. However, in practice one is interested primarily in the 
lowest-lying levels; and we find significant evidence that the Bethe Ansatz solution yields 
the ground state energy with just k = 1 (i.e., M = |iV) for iV even, and with k = (i.e., 
M = i(iV — 1)) for iV odd. We also investigate the special case 

a_ = -a+ , p- = -P+ , 6+ = 9_ = , iV = odd (1.9) 

considered in j^j, and find evidence that the Bethe Ansatz solution gives the complete set of 
eigenvalues with k — 0. 

The outline of this article is as follows. In Section we consider the Bethe Ansatz 
solution for generic values of both bulk and boundary parameters. We describe an ingenious 
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procedure, which was pioneered by Barry McCoy and his collaborators, for addressing the 
problem of completeness; and we present the results we have obtained using this procedure. 
In Section 01 we consider the special case (jl.9|) . We conclude with a summary of our main 
results in Section HJ 

2 Generic case 

We begin by addressing the following question: given a value of N, what value of k is 
needed to obtain from the Bethe Ansatz solution (jl.2|) - (jl.8|) the complete set of energy 
eigenvalues? To clarify the meaning of this question, let us make the elementary observation 
that the number of Bethe roots must satisfy M > 0; hence, (II .8|) implies that k must satisfy 
k > 1 — N. The minimum value k = 1 — N corresponds to zero Bethe roots, and therefore, 
to only one eigenvalue. (See Eq. ()1.3JI .) Since a chain with N spins has 2 N eigenvalues 
(which are distinct for generic values of parameters), this minimum value of k can give a 
complete set of eigenvalues only for N = 0. 3 We wish to determine, for higher values of N, 
the value (s) of k needed to obtain from the Bethe Ansatz solution the complete set of 2 N 
eigenvalues. 

Such questions of completeness are notoriously difficult to address, even numerically. 
Indeed, since there is no known systematic way of solving Bethe Ansatz equations, it is 
not possible to decide unequivocally when one has found all the solutions of those equations. 
Fortunately, there does exist a systematic method, exploited by McCoy and his collaborators 
(see, e.g. jHl [TDJ ) , of determining the Bethe roots corresponding to a given eigenvalue. 
Since, for small values of N, the eigenvalues can be computed by direct diagonalization, 
this method can be used to determine whether the Bethe Ansatz solution reproduces all the 
known eigenvalues. (Since in this approach one does not actually "solve" the Bethe Ansatz 
equations, the possibility remains open that those equations may admit additional solutions 
which do not correspond to actual eigenvalues. We shall return to this point later in Section 

E3) 

This method, to which we refer as 'McCoy's method', actually makes use of the full 
transfer matrix of the model, rather than the Hamiltonian. Hence, we now briefly review its 
construction. 

3 Although the Hamiltonian makes sense only for N > 2, the transfer matrix (which is described 

below) is well-defined even for N = 0. 



3 



2.1 Transfer matrix 



The transfer matrix for the open chain (Jl.lj) is constructed according to Sklyanin's recipe j2] 
from the R matrix 



R{u) 



( sinh(w + 77) \ 

sinh it sinhr; 

sinh77 sinh it 

\ sinh(u + 17) ) 



(2.1) 



and the 2x2 nondiagonal matrices K T (u) whose components are given by j3*j ITT] 

K-ii ( u ) = 2 (sinh a_ cosh j3_ cosh u + cosh a_ sinh /3_ sinh it) 

i^22 (it) = 2 (sinh a_ cosh /3_ cosh u — cosh a__ sinh /3_ sinh it) 

K^ 2 {u) = e e ~sinh2it, K 2X (it) = e _e ~ sinh 2it , 

and 



(2.2) 



-2 (sinh a + cosh (3+ cosh(it + 77) — cosh a + sinh /3 + sinh (it + 17)) 
-2 (sinh a + cosh /3+ cosh(ix + 17) + cosh a + sinh (3 + sinh(it + 77)) 
-e 9+ sinh2(w + 77) , K^(u) = -e~ 0+ sinh 2(u + 17) . (2.3) 



The matrices X =F (it) are equal to those appearing in |Hj divided by the factors k t , respectively. 
This leads to the rescaling of the transfer matrix already mentioned in Footnote 2. 

The transfer matrix t(u) is given by [2j 



t(u) = tr K+(u) T (u) K (u) f (u) 
where the monodromy matrices are given by 



T (it) = RonM ■ ■ ■ Roi(u) 



T (it) = RqAu) ■ ■■R QN (u) 



(2.4) 



(2.5) 



and tro denotes trace over the "auxiliary space" 0. The transfer matrix has the important 
commutativity property 



[t(u),t(v)} = 0. 
and it "contains" the Hamiltonian (jl.lj) . 



u=0 



+ C2I. 



(2.6) 



(2.7) 



where 

Cl 



16 sinh a_ cosh (3- sinh a + cosh j3 + smh 2N 1 77 cosh rj ' 



sinh 2 1] + N cosh 2 77 ^ 
2 cosh 77 

and I is the identity matrix. According to the Bethe Ansatz solution j^j, the eigenvalues 
A(w) of the transfer matrix are given by 

m = h(u) 9^=A + h( ^- n) 9^A } (2 . 9) 

Q{u) Q(u) 

where h(u) and Q(u) are given by (jl.5j) and (jl.fi j) . respectively. The formula (jl.Hj) for the 
energy eigenvalues follows directly from (j2.7j) - (j2.9j) . 



2.2 McCoy's method 

To implement McCoy's method, it is more convenient to work with the spectral parameter 
x = e u and the anisotropy parameter q = e v . We denote by t(x) the transfer matrix expressed 
in terms of x, and similarly for other quantities. 

McCoy's method consists of four steps: 4 

1. Fixing an arbitrary value x$ of the spectral parameter, compute numerically the eigen- 
vectors I A) of the transfer matrix t(xo). Due to the commutativity property of the 
transfer matrix, the eigenvectors do not depend on the spectral parameter. 

2. Determine the eigenvalues A(x) by acting with t(x) on the eigenvectors |A). Due 
to the commutativity property of the transfer matrix, these eigenvalues are Laurent 
polynomials in x. 

3. Set Q(x) = J2h=o^k (x 2k + (xq)~ 2k ) (see (jl.fi j) ). and determine the coefficients {bk} 
from the relation (123, i.e., A(x)Q(x) = h(x)Q(2) + h(\)Q(xq). 

4. Factor the polynomials Q(x), whose zeros Xj are the sought-after Bethe roots. 

Below we present the results that we have obtained using this method. We discuss 
separately the "massless" regime (77 is purely imaginary) and the "massive" regime (77 is 
purely real). 

4 We are grateful to B. McCoy for explaining this procedure to us. 
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2.3 Massless regime 



For the case that the bulk anisotropy parameter 77 is purely imaginary, the transfer matrix 
is generally not Hermitian. Hence, in principle, it may have fewer than 2 N eigenvectors. 
Nevertheless, if the boundary parameters are suitably restricted, the transfer matrix can be 
shown to be a normal matrix, i.e., 

[t(u), t(«)t] =0, (2.10) 

which implies that it is unitarily diagonalizable. Indeed, treating the spectral parameter u as 
real, it is easy to see that the R matrix (j2.1j) satisfies R(uy = —R(—u). Let us now restrict 
the boundary parameters so that 

a T , 9 T = purely imaginary , (3 T = purely real . (2-H) 

The K matrices ()2.2j) . (j2.3J) then obey similar relations K T {u) J[ = —K T (—u). It follows that 
the transfer matrix obeys the simple relation 

t( U y = t(-u). (2.12) 

Combining this result with the commutativity relation (j2.fi)) immediately yields the desired 
result (|2.1())1 . Moreover, the conditions (|2.11)1 for the boundary parameters imply that the 
Hamiltonian ()1.1|) is Hermitian. 

In the numerical work which we present below, the values of the bulk and boundary 
parameters are chosen as follows: 

7] = 0.3i, a + = 0.75i, (5 + = -0.5 , 6 + = 0.8i , 
a_ = 0.25i+(fc-l)(0.3i), ^- = 0.5, 0_ = O.li. (2.13) 

This set of values satisfies the constraint p. 2)1 for any value of k, as well as (j2.11j) . 

Tabled shows, for values of N ranging from to 4, all the 2 N energy eigenvalues and the 
corresponding shifted Bethe roots. Our main observation is that, for each value of N, the 
corresponding value of k is equal to N + 1. (We obtained similar results for up to iV = 7, 
but we do not present the data here.) For k > N + 1, the Bethe Ansatz also yields all 2 N 
energy eigenvalues. But for 1 — N < k < N + 1, the Bethe Ansatz does not yield all 2 N 
energy eigenvalues. 5 In other words, the minimum value of k for which the Bethe Ansatz 
reproduces all the eigenvalues is k = N + l. We have observed this numerically for the choice 

5 For the "missing" eigenvalues (i.e., those eigenvalues of the transfer matrix which are not given by the 
Bethe Ansatz solution), step 3 of McCoy's method fails: for such eigenvalues A(x), there are no appropriate 
polynomials Q(x) which satisfy A(x)Q(x) — h(x)Q(^) + h(^-)Q(xq). 
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of parameters ()2.13|) with values of N up to 7, and we conjecture that it is true for generic 
values of the boundary parameters for all N. 

We further conjecture that for k > N + 1, the Bethe Ansatz equations (jUJ) admit extra- 
neous solutions, which do not correspond to eigenvalues of the transfer matrix. We have 
verified this numerically for N = with k = 3. Indeed, with the choice of boundary pa- 
rameters (j2.13|) . we find for this case not only the Bethe root u = 0.690849 + 1.57082 which 
corresponds to the (single) transfer matrix eigenvalue, but also an additional solution of 
the Bethe Ansatz equation u = 1.4208i which does not correspond to this eigenvalue. As 
emphasized in the beginning of Section |2l it is difficult to hunt for solutions of Bethe Ansatz 
equations, especially for higher values of N and M. 

Assuming that these two conjectures are correct, it follows that the Bethe Ansatz yields 
all 2 N eigenvalues and no extraneous solutions for precisely k = N + 1. While it is gratifying 
to obtain all the eigenvalues, this result is also disappointing, since the constraint (jl.2j) then 
implies that the imaginary parts of the boundary parameters should grow linearly with N. 

2.3.1 Ground state 

Although a high value of k is required to obtain all the energy levels (namely, k — N + 1), 
we find that the ground-state energy can be obtained with a much lower value of k. Indeed, 
using the parameter values (j2.13J) . we performed a search for the minimum value of k (for a 
given value of N) where the Bethe Ansatz reproduces the ground-state energy, up to N = 7. 
Our results are summarized in Table which gives in addition to the value of k also the 
ground-state energy and the corresponding Bethe roots. Our main observation here is that 
k = for N odd, and k = 1 for N even. We conjecture that this result is true for generic 
values of the boundary parameters for all N. If correct, then the Bethe Ansatz is practical 
for investigating the ground state in the thermodynamic limit. 

We also observe from Tabled that the shifted Bethe roots are real for the ground state, 
as is also the case for the closed XXZ chain with periodic boundary conditions. (For higher 
values of k, the shifted Bethe roots for the ground state are either real or have imaginary 
parts «7r/2, as can be seen from Tabled) 

Finally, we remark that our numerical results suggest that the Bethe Ansatz correctly 
yields 2 N ~ 1 eigenvalues for k = (N odd), and 2 N ~ 1 + |(jw 2 ) eigenvalues for k — 1 (N 
even). 6 

6 In formulating the latter conjecture, which we have checked up to N = 8, a useful reference was 
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2.4 Massive regime 



For the case that r\ is purely real, we choose 0q= = and the remaining boundary parameters 
to be real, thereby making the transfer matrix manifestly Hermitian. In particular, for the 
numerical work presented below, we take the values 

v = 0.3, a+ = 0.75, (3 + = -1.2, 9 + = , 
a_ = 0.25 + (k -1) (0.3), /3_ = 0.5 , 9_ = , (2.14) 

which satisfy the constraint (jl.2|) for any value of k. 

Our results for the massive regime are very similar to those for the massless regime. 
Indeed, consider Table El which shows all the 2 N energy eigenvalues and the corresponding 
Bethe roots for values of iV ranging from to 4. As in the massless case, the minimum value 
of k for which the Bethe Ansatz reproduces all the eigenvalues is k = N + 1. (We obtained 
similar results for N = 5. For larger values of N, roundoff errors become significant.) 
Moreover, as shown in Table IH the Bethe Ansatz reproduces the ground state energy for 
k = for N odd, and k = 1 for N even. The corresponding shifted Bethe roots are purely 
imaginary. 



3 Special case 

We now turn to the special case (jl.9|) . which was first considered in (SJ- in this case, the 
boundary terms of the Hamiltonian (jl.lj) reduce to 

- sinh 1] [cotha_ tanh/3_ (af — a z N ) + csch ct_ sech/3_ (cr^ — a x N )} . (3.1) 

We first argue that for this case all the energy eigenvalues are 2-fold degenerate. Indeed, it 
is easy to see that the Hamiltonian commutes with the operator U defined by 7 

U = C P, (3.2) 

where C is the "charge conjugation" operator 

N 

C=l[al, (3.3) 

n=l 

7 A similar symmetry operator was invoked in to argue that an open chain with diagonal boundary 
terms has a two-fold degenerate spectrum for N odd. There the argument is simpler, since in that case the 
Hamiltonian also commutes with 5* z , while U and S z anticommute. 
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which satisfies = C and C 2 = 1; and P is the "parity" operator ^3], which satisfies 

PalP = 4 +1 _ n , (3.4) 

as well as P^ = P and P 2 = 1. It follows that also U is Hermitian and squares to 1. Hence, 
U has eigenvalues ±1. For N odd, U has an equal number of +1 and —1 eigenvalues. 8 It 
follows that all energy eigenvalues are 2-fold degenerate. In fact, since U commutes with the 
full transfer matrix t(u), all the eigenvalues A(u) are 2-fold degenerate. 

Since for this case there are generally only 2 N ~ l distinct eigenvalues, one expects that 
all of these eigenvalues can be reproduced by the Bethe Ansatz with a value of k < N + 1. 
Indeed, as shown in Table El we find significant evidence which supports the conjecture that 
the complete set of 2 iV_1 eigenvalues is obtained for k = (i.e., M = h(N — 1)). (We 
obtained similar results for N = 7.) 



4 Conclusion 

Within the range of parameters which we have explored (as detailed in Sections 12.31 and 
12.4)1 . we have found significant numerical evidence for the following conjectures regarding 
the Bethe Ansatz solution JOJ) - dTTHJ) of the model (TTTJ) : 

• For generic values of the bulk and boundary parameters satisfying (jl.2j) . the solution 
yields the complete set of 2^ eigenvalues for k = N + 1 (i.e., M = N). 

• The solution yields the ground state energy for k = 1 (i.e., M = ±N) when N is even, 
and for k = (i.e., M — ~(N — 1)) when iV is odd. In the massless regime, the shifted 
Bethe roots corresponding to these states are real. 

• In the special case ()1.9)) where the spectrum is 2-fold degenerate, the Bethe Ansatz 
solution yields the complete set of 2 N ~ l eigenvalues for k = (i.e., M = |(JV — 1)). 

These results suggest that the Bethe Ansatz solution is both valid and practical for 
investigating the ground state (and presumably, also low-lying excited states) of the model 

8 To prove this, it suffices to show that the trace of U is zero. For N odd, the parity operator leaves the 
"middle" spin at site ^(N + 1) invariant. Hence, 

tr U = tr 12 ... N U = tr i (N+1) (^| (JV+1) ) tr' (p J] o*) = , 
since the Pauli matrix a v is traceless. (Here tr' denotes trace over all spaces n ^ \{N + 1).) 
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(jl.lj) in the thermodynamic limit. In particular, these results provide justification for the 
computations in ^H] of the thermodynamic limit for the special case (jl.9j) . and clear the 
way for analogous computations in the general case. We stress, however, that this model 
has many parameters, other ranges of which remain to be explored. 
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Note added: 

There is a counterpart of the special case (jl.9|) for even values of N, namely, 

«_ = — a + + 7] , P- = —(3+ , 6+ = 0- = , N = even , 

and hence k — 1. For this case the spectrum also has degeneracies. For even values of N up 
to N = 6, we find that the Bethe Ansatz solution with M = |iV gives the complete set of 
2-W-i _|_ ±(^ 2 ) distinct eigenvalues. 
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ADDENDUM to "Completeness of the Bethe Ansatz solution of the open 
XXZ chain with nondiagonal boundary terms" 

In [1] (to which we refer hereafter by I), we find significant numerical evidence that for 
the open XXZ quantum spin chain Hamiltonian (II. 1) with bulk and boundary parameters 
satisfying the constraint (11.2) 

a_ + (3. + a + + f3 + = ±(6- - 9+) + rjk , (1) 

the Bethe Ansatz solution (11.3) - (11.8), (12.9) gives the complete set of 2 N eigenvalues for 
k = N + 1; and this solution gives only some (but not all) of the eigenvalues for 1 — iV < k < 
N + 1. Here we conjecture that a simple generalization of this Bethe Ansatz solution gives 
the complete set of eigenvalues for all values of k in the interval — (N + 1) < k < N + 1. 

Indeed, consider the following two expressions, distinguished by ±, for the eigenvalues of 
the transfer matrix (12.4): 

A±(u) = h±{u) +h±{-u-ri) , 2 

Q±{u) Q±{u) 

where h±(u) are given by 



h±(u) = — sinh (u + r]) — 



2N , . sinh(2w + 277) 



sinh(2w + rj) 

x 4 sinh(u ± a J) cosh(w ± (3-) sinh(w ± a+) cosh(w ± f3 + ) ; (3) 



and Q±(u) are given by 



where 



M± 

Q±( u ) — Y\_ s i n h( u — u f) sinh(u + uf + if) , (4) 



M ± = 1 -{N - 1 ± k) . (5) 

The corresponding Bethe Ansatz equations are 

h±{uf) Q±(uf + V ) . 
Tl — T 7 = -7TT^E r> j = l,...,M±. (6) 

h±{-u j - 7]) Q±{u j - 1]) 

The Bethe Ansatz solution considered earlier [2,1] corresponds to the above expressions with 
the plus (+) sign. 

We conjecture that for a given set of bulk and boundary parameters satisfying (JJJ with 
\k\ < N + 1 and k odd (even) for N even (odd), the eigenvalues A + (u) and A_(u) together 



12 



give the complete set of eigenvalues of the transfer matrix. That is, for k = N + 1, all the 
eigenvalues are given by A + (u) (as found in I); for k = —(N + 1), all the eigenvalues are 
given by A_(u); and for —(N + 1) < k < N + 1, both A + (u) and A_(w) are needed to obtain 
the complete set of eigenvalues. 

This conjecture is supported by significant numerical evidence, obtained (as in I) using 
'McCoy's method'. For instance, for N = 6 and k — 1, we find that 42 eigenvalues are given 
by A + (u) and 22 eigenvalues are given by A_(u); together, they give the complete set of 
2 6 = 64 eigenvalues. 
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N 


k 


M 




shifted Bethe roots Uj 





1 





0.259617 


- 


1 


2 


1 


-0.455662 


0.278824 








0.455662 


0.702627 + u 


2 


3 


2 


-1.48624 


0.124133, 0.691313 + u 








0.13202 


0.334031 ± 0.186753i 








0.483745 


0.467053, 0.678121 + u 








0.870471 


0.55974 + M , 0.974697 + u 


3 


4 


3 


-2.20345 


0.0866988, 0.373369,0.655889 + u 








-1.75971 


0.0794522, 0.542856 + m , 0.939102 + m 








-0.0127446 


0.240758 ± 0.139792i , 0.659713 + u 








-0.00644228 


0.21088, 0.541397 + m , 0.935445 + m 








0.671117 


0.39724, 0.362185 ± 0.341578* 








0.878343 


0.638761 + M , 0.511044 ± 0.205038i 








1.09747 


0.615979 , 0.530429 + u , 0.906019 + u 








1.33542 


0.468162 + M , 0.719202 + m , 1.16451 + u 


4 


5 


4 


-3.20655 


0.0656969 , 0.173673 , 0.505631 + m , 0.869254 + m 








-2.15365 


0.0641635, 0.601829 + m , 0.390452 ± 0.193885i 








-1.895 


0.0609306 , 0.526013 , 0.497497 + m , 0.847927 + m 








-1.6139 


0.0587044 , 0.432076 + u , 0.672212 + u , 1.10249 + u 








-0.66516 


0.163048, 0.365317 ± 0.193132z, 0.602513 + u 








-0.549331 


0.145871 , 0.515342 , 0.497186 + u , 0.847232 + u 








-0.345969 


0.136943 , 0.431537 + u , 0.671367 + u , 1.101 + u 








0.259969 


0.506611 + M , 0.87182 + m , 0.160443 ± 0.150348i 








0.842045 


0.287219, 0.297174 ± 0.310653? , 0.609213 + m 








0.891361 


0.497959 + M , 0.849881 + m , 0.370167 ± 0.118987i 








0.948984 


0.285065 , 0.429407 + w , 0.668011 + m , 1.09495 + w 








1.1934 


0.429338 ± 0.160247i, 0.373211 ± 0.475013i 








1.33834 


0.56912, 0.582163 + u , 0.533392 ± 0.368923i 








1.48584 


0.48142 + M , 0.807965 + m , 0.649061 ± 0.215203i 








1.64707 


0.733488 , 0.41686 + m , 0.647858 + u , 1.05404 + m 








1.82255 


0.367461 + M , 0.570419 + m , 0.839668 + m , 1.29177 + w 



Table 1: Complete set of 2 energy levels and Bethe roots in the massless regime, using pa- 
rameter values (|2.13J) . We use the shorthand notation «o = ivr/2. Without loss of generality, 
we restrict the shifted Bethe roots so that 3?e Uj > and — | < 3m v,j < f . 
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N 


k 


M 


ground-state energy E 


shifted Bethe roots Uj 


1 








-2.79413 




2 


1 


1 


-1.57715 


0.0944455 


3 





1 


-4.58216 


0.0973252 


4 


1 


2 


-3.3477 


0.0559452, 0.139137 


5 





2 


-6.35177 


0.0572972, 0.14122 


6 


1 


3 


-5.10816 


0.0402978, 0.0893334, 0.168789 


7 





3 


-8.11181 


0.0410605, 0.0906602, 0.170368 



Table 2: Ground-state energy and Bethe roots in the massless regime, using parameter values 
(EES)- 
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AT 

iv 


7 

fc 


?\ /f 
M 


Z7 1 


shifted Bethe roots Uj 





1 


u 


a 0001 r 

(J. 28216 




1 


2 


l 


-0. 478182 


0.274302? 








a a to i on 

U. 478182 


"1 OOOOT i 

1.82287 + u 


2 




2 


-1.52836 


a i o o o a/ 1 ' i oo/io i 

0.123806z, 1.83483 + u 








0.128052 


r\ i o /ion/ 1 i a o oto/ 1 r " 

0.184296 ± 0.327865^ 








0. 496614 


A A £?AT£?0 * 1 O /I A/ 1 1 i 

0.460762?, 1.84961 + Mo 








a aao/^ao 

0.903692 


1.44453 + Mo , 2.2506 + Mo 


3 


4 


3 


O OTO O 

-2.2733 


a r\o r* r* >^ try ' a o/^a/' 1 a a * i ot/Ioa i 

0.0866767? , 0.369644? , 1.87489 + m 








-1.80933 


a a 7A /l A o A ' 1 /irorr i o o T a a 1 i 

0.0794024?, 1.45356 + Mo, 2.27091 + Mo 








-0.00928535 


0.140292 ± 0.238892? , 1.86994 + Mo 








r\ A A /I Z" 1 O O O/" 1 

-0.00463236 


a o i a r o o ' i a r a r c\ i oo to o o i 

0.210583z, 1.45452 + uq , 2.27283 + u 








0.684253 


A O O TO O O ' A O O TO A T 1 A OT T1 O O ' 

0.387332z, 0.337297 ± 0.355183? 








0.903368 


AOAO/1T/1 1 A r ATAO O ' A O A TO O i 

0.202474 ± 0.507083?, 1.89723 + m 








1 1 O "7 1 ' ) 

1.12743 


A C* A0 1 A A ' 1 /l/^O/ir i OO O TT/ 1 i 

0.608104?, 1.46245 + Mo , 2.28776 + Mo 








1.38149 


1 OTO /I O i "1 TO /1/1T i OTAAOO i 

1.36248 + Mo , 1.78447 + Mo , 2.59988 + Mo 


A 

4 


5 


A 

4 


-3.30164 


A A/^CTAOT^' A 1 TOTO/'V 1 /1TOOA I n O 91 /1AO I 

0.0657087?, 0.173726? , 1.47329+ Mo, 2.31403 + Mo 








-2.21969 


A (\0 A A O A A ' A 1 AO A C 1 1 A O O O A "1 /I ' 1 C\ A ^ A H i 

0.0641814?, 0.192461 ± 0.388014? , 1.94346 + m 








-1.95389 


A C\CC\C\CCC\A. A CO/1/1AA^ 1 A TAAT I ~ O OOCO^ I 

0.0609669?, 0.524409? , 1.47997 + Mo, 2.32536 + Mo 








1 / ' / ■ O *" O 

-1.66252 


A ATOTOA/" 1 ' 1 O/^TAI i 1 TAAA"1 i O 1 ' O / " *" ' ) i 

0.0587306z, 1.36591 + u , 1.79991 + u , 2.62653 + u 








-0.68481 


A 1 A 1 T /I Z" 1 1 AO/ 1 OTAO ' A 1 Z" 1 O O O O ' 1 A /I OO A i 

0.191746 ± 0.362798?, 0.163232? , 1.94239 + Mo 








r\ r st ta 

-0.56579 


A 1 /I Z" 1 A A T ' A f 1 O A O T ' "1 A OAO i O O O f TZ" 1 i 

0.146007z, 0.513987? , 1.4802 + m , 2.32576 + m 








-0.355209 


a i otao/1^' i orcnn i „ ioaaoo i ~ o z^otai i 

0.137034z, 1.36599 + u , 1.80022 + Mo, 2.62701 + Mo 








0.270674 


n 1 rnon i a 1 rnn o 1 a a to r'Oi o o i o z" 1 a i 

0.150351 ± 0.160618?, 1.47252+ m , 2.31269 + m 








a o to r /i 

0.87254 


A01AOOT 1 AOATOOO' AOOrOO/1' 1 AO "1 TZ" 1 i 

0.310225 ± 0.295283? , 0.285824? , 1.93176 + Mo 








0.920906 


A11AATO 1 A O T A TO ' 1 A T A O 1 i O OO A A O i 

0.119052 ± 0.37073? , 1.47931 + m , 2.32448 + m 








a atootc 

0.978875 


n noc^i c i o z? z? o o i ~ i oni ho i ~ o z? o o a 1 i 

0.285415? , 1.36633 + Mo , 1.80148 + Mo , 2.62891 + Mo 








l.zzzlo 


U.iooifo ± U.414(l? , U.4t)oyo ± U.oo40oOZ 








1.38081 


0.366909 ± 0.533061H , 0.565907?, 1.97405+ m 








1.52904 


0.212615 ± 0.654327? , 1.49436 + u , 2.34875 + m 








1.68784 


0.724014* , 1.36869 + u , 1.80952 + m , 2.64055 + m 








1.8807 


1.35068 + u , 1.66925 + u , 2.10024 + u , 2.91917 + m 



Table 3: Complete set of 2 energy levels and Bethe roots in the massive regime, using pa- 
rameter values ()2.14|) . We use the shorthand notation uq = in/2. Without loss of generality, 
we restrict the shifted Bethe roots so that 9fte iij > and — f < Qro Uj < |; or 5Re Uj = 
and < iij < |. 
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N 


k 


M 


ground-state energy E 


shifted Bethe roots Uj 


1 





o 


-2.86459 




2 


1 


l 


-1.61648 


0.0941058i 


3 





l 


-4.71323 


0.0969713i 


4 


1 


2 


-3.4432 


0.0558492i, 0.138676i 


5 





2 


-6.53887 


0.0571973i, 0.140752? 



Table 4: Ground-state energy and Bethe roots in the massive regime, using parameter values 
(EH- 



AT 

IV 


h 
K 


IVl 


TP 
Hi 


shifted Bethe roots iij 


1 


U 


U 


U 




o 

O 


U 


1 


-Z.U606I 


U.U8llz8 I 








n ononQO 
-U.zUz9o8 


C\ 00QQ70 








U.98ol6 I 


i.i/ /y + % 








1 OSQQS 


U.00 / \j06l 


c 



U 


o 
L 


Q OQO/1 Q 

-o.yoz4o 


n nci 0070 n 101 Qn/l 

U.Uoloo/o, U.l/loU4 








-z.ooouo 


n n^ni ^91 n 9^qq&7 








-1.61887 


0.0470911, 1.18204 + u Q 








-1.53627 


0.114268, 0.255101 








-1.3471 


0.0452599, 0.570638i 








-0.687493 


0.10434, 1.18027 + u 








-0.474743 


0.0991988, 0.56857i 








0.268219 


0.119246 ± 0.149991i 








0.486895 


0.19328, 1.17489 + u 








0.642915 


0.17875, 0.564459i 








0.954064 


0.278324 ± 0.155511z 








1.49874 


0.407709, 1.14897 + u 








1.65895 


0.347419 , 0.557976* 








1.96938 


0.850516 + w , 1-62572 + u 








2.28162 


0.551706i, 1.24328 + u 








2.3948 


0.550319i, 0.950981i 



Table 5: Complete set of 2 Ar_1 energy levels and Bethe roots for the special case (jl. 9)1 . with 
r] = 0.3i, a_ = — a + = OAi, (3- = —(3+ = 0.7, 9 + = 9- = 0. We use the shorthand notation 
Uq = in/2. Without loss of generality, we restrict the shifted Bethe roots so that 3?e Uj > 
and — I < 3m itj < |; or 3?e itj = and < 3m itj < |. 
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